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Abstract 

We present two classes of random walks restricted to the quarter plane whose gen- 
erating function is not holonomic. The non-holonomy is established using the iterated 
kernel method, a recent variant of the kernel method. This adds evidence to a recent 
conjecture on combinatorial properties of walks with holonomic generating functions. 
The method also yields an asymptotic expression for the number of walks of length n. 
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1 Introduction 

Previous studies of random walks on regular lattices had led some to conjecture that random 
walks in the quarter plane should all have holonomic generating functions. This is not 
unreasonable, given that walks confined to the half plane have algebraic generating functions 
(PQ proves this for directed paths, but their results extend to random walks). This was, 
however, disproved by Bousquet-Melou and Petkovsec with their proof that knight's walks 
confined to the quarter plane are not holonomic pE]. 

In this paper we give two new examples of random walks whose generating functions are 
not holonomic. We use a technique developed in a recent study of self-avoiding walks in 
wedges [H] called the iterated kernel method. This is an adaptation of the kernel method (2j 
[TO] , in which we express a generating function in terms of iterates of a kernel solution. We 
can show that this generating function has an infinite number of poles, which is incompatible 
with holonomy. 

The aim of this work is two-fold. The first goal is to illustrate a potentially general 
technique to prove non-holonomy. Secondly, this work was completed in the context of a 
general generating function classification of all nearest neighbour lattice walks [9], and lends 
evidence to a general conjecture on the connection between different symmetries of the step 
sets and the analytic nature of their generating functions. 
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Figure 1: Sample walks with steps from S = {NW, NE, SE} and T = {NW, N, SE} 



1.1 Walks and their generating functions 

The objects under consideration are walks on the integer lattice restricted to the quar- 
ter plane with steps taken from S = {(— 1, 1), (1, 1), (1, — 1)} in the first case, and T = 
{( — 1, 1), (0, 1), (1, —1)} in the second case. We use the compass notation and label the di- 
rections of these two sets as S = {NW, NE, SE} and T = {NW, N, SE}. Two sample walks 
are given in Figure 11.11 

We associate to these two steps sets two power series each: W(t) a counting (ordinary, 
univariate) generating function and Q(x,y;t) a complete (multivariate) generating function. 
The series W(t), is the ordinary generating function for the number of walks. It is a formal 
power series where the coefficient of t n is the number of walks of length n. The complete 
generating function Q(x, y; t) encodes more information. The coefficient of x l yH n in Q(x, y; t) 
is the number of walks of length n ending at the point Remark that the specialization 

x — y — 1 in the complete generating function is precisely the counting series, i.e. Q(l, 1; t) = 
W(t). If the choice of step set is not clear, we add an index. 

In part, our interest in the complete generating function stems from the fact it satisfies 
a very useful functional equation which we derive using the recursive definition of a walk: a 
walk of length n is a walk of length n — 1 plus a step. The quarter plane condition asserts 
itself by restricting our choice of step should the smaller walk end on an axis. 

The step set S leads to the following equation: 

Q(x, y;t) = l+t(xy + - + ^] Q(x, y; t) - t-Q(x, 0; t) - t^Q(0, y; t), (1.1) 
\ y xj y x 

and set T defines a similar equation: 

Q(x, y;t) = l+t(y+- + ^] Q(x, y; t) - t-Q(x, 0; t) - ^Q(0, y; t). (1.2) 
V y xj y x 

1.2 Properties of holonomic functions 

We are interested in understanding the analytic nature of the generating functions. This 
nature gives a basic first classification of structures, and also some general properties such as 
general form of the asymptotic growth of the coefficients. See for example Bousquet-Melou's 
recent work on classifying combinatorial families with rational and algebraic generating func- 
tions [3j. We are interested in generating functions which are holonomic, also known as 
D-finite. Let x = Xi, x 2 , ■ ■ ■ , x n . 
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Definition 1.1 holonomic function. A multivariate function G(x) is holonomic if 
the vector space generated by the partial derivatives of G (and their iterates), over rational 
functions of x is finite dimensional. This is equivalent to the existence of n partial differential 
equations of the form 



for i satisfying 1 < i < n, and where the pj^ are all polynomials in x. 

The basic feature of these functions that we use is the fact that they have a finite number 
of singularities. For example, in the univariate case, these are poles which appear are zeroes 
of p dl ,x- 

1.3 Main results 

We show that both the complete and counting generating functions of the walks with steps 
from S = {NW, NE, SE} and T = {NW, N,SE} are not holonomic using this singularity 
property of holonomic functions. In fact, it suffices to show that their counting generating 
function W(t) is not holonomic, as holonomic functions are closed under algebraic substitu- 
tion [8j [13], and we have W(t) = Q(l, l;i). We show that W(t) has an infinite number of 
poles and thus is not holonomic [13]. This will give the first main theorem that we present. 

Theorem 1.1. Neither the complete generating function Qs(x,y;t) , nor the counting gener- 
ating function Ws(t) = QsO-i 1; t) of nearest-neighbour walks in the first quadrant with steps 
from {NE, SE, NW} are holonomic functions with respect to their variable sets. 

To prove this, we use the iterated kernel method, as described by Janse van Rensburg et 
al. in [IT] in counting partially directed walks confined to a wedge. 

The basic idea is similar to other variants of the kernel method in that we write our 
fundamental equation in the kernel format, determine particular values of x and y which fix 
the kernel, and then generate new relations. For these walks, the kernel iterates form an 
infinite group whose elements each introduce new poles into the generating function. The 
main difficulty in proving that these generating functions are not holonomic is proving that 
the singularities do not cancel. 

As step set T is not symmetric in the line x = y, the relations obtained are more complex 
and while the argument is essentially the same, the details of the proof of its non-holonomy 
is more complicated, as we see in Section [31 This will give the second main theorem. 

Theorem 1.2. Neither the complete generating function Qr{x,y;t) , nor the counting gen- 
erating function Wr{t) = <3t(1,1;<0 of nearest-neighbour walks in the first quadrant with 
steps from {NW, N, SE} are holonomic functions with respect to their variable sets. 
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2 The non-holonomy of Qs{x, y] t) 

We use the iterated kernel method to find a closed form expression for the generating func- 
tion Qs(x,y;t). We then analyse this expression to prove this function, and the counting 
generating function Ws(t) are not holonomic. 

2.1 Defining the iterates of the kernel solutions 

To begin, consider the kernel version of Equation (II. lft : 

(xy - tx 2 y 2 - tx 2 - ty 2 ) Q(x, y) = xy - tx 2 Q(x, 0) - ty 2 Q(y, 0). (2.1) 

Here, for brevity we write Q(x, y \ t) as Q(x, y), and have used the x y symmetry to rewrite 
Q(0,y) asQM). 

Since the kernel K(x, y) = xy — tx 2 y 2 — tx 2 — ty 2 is a quadratic, it has two solutions as 
a function of y 

Y ±1 (x) = 2t{i X +x2) (l =F Vl- 4*3(1+**)) . (2.2) 
Note that as formal power series in t, these roots are: 

Y +1 (x) = xt + 0(t 3 ) 

Y^{x) = -^--xt + 0{e). 
(1 + x z ) t 



Further, Y + i(x) is a power series in t with polynomial coefficients in x. We will need to form 
compositions of these roots. Write Y n (x) = (Yio) n (x), and likewise for Y^ n . 

Lemma 2.1. The kernel roots obey the relation: 

r +1 (y_i(x)) = x and = x. (2.3) 

and the set {Y n \n e Z} forms a group, under the operation Y n (Y m (x)) = Y n o Y m = Y n+m , 
with identity Yq — x. Furthermore, they also obey the relation: 

111 . , 

+ T? 7Z\ = ZT> ( 2 - 4 ) 



Yx{x) Y-xix) tx 
which extends (by substituting x = Y n -x(x)) to: 

1 1 



)],(.r) *y n _! Y n . 2 (x)- (2 ' 5) 
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Proof. To prove Equation (12. 3p . consider the four compositions of the kernel roots Y±(Y±(x)) 
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and denote them by {C{\. We form the product TT(C — Ci)j f° r formal parameter C, and 



i=i 



one may verify: 

(x - C) 2 ((t 2 + x 2 )C - x(l - 2t 2 )C + t 2 x 2 ) = 0. 

Hence two of the compositions are the identity. 

There are L) possibilities, and we can exclude 4 of these by consistency arguments - 
i.e. they will imply that Y + = K_ which is not true. This leaves us to check that either 
Y + o Y_ = Y_ o Y + = x or Y + o Y + = Y_ o Y_ = x. A manual calculation then shows that 
Y_ o Y + = Y + o Y_ = x. 

Equation (12.41) follows since Y±i are roots of a quadratic — in particular Y + {Y_i and 
Y + i + y_i are rational functions of x and t. □ 

2.2 An expression for Q(x, 0) in terms of the iterates 

We are now able to find a formal expression for the generating function. 

Theorem 2.2. The generating function Q(x,0) is given by 

Q(x,0) = 4t V(-l) n r n (x)Wx). (2-6) 

and so Q(l, 1) is given by 

l-2 f Q(l,0) = l-2 E „, (-l)-Y„(l)^, +1 (l) 
^ V ' ; l-3t 1-3* K J 

Proof. By construction we have K(x,Yi(x)) = 0, thus substituting y = Y\(x) into Eq. (12.11) 
gives (after a little tidying) 



^o) = ^-^«),o). 

a; r 



Now we substitute x = F n (x) into this equation to obtain 

^•°»Kw)^(w) 2gM0) ' 

Using this expression for Q(Y n (x),0) for increasing n, we can iteratively generate a new 
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expression for Q(x,0): 



x t x 2 \Y 1 (x)t Y^x) 2 

M X )Y 1(X) u x) y 2 ( x ) + 



xH xH \ x 



N-l f ^ r ^__\\ 2 

Q(Y N (x),0). 



N I Yn{%) 



1 / 

= ^E(- 1 ) n ^)^^) + (- 1 ) JV ,. 

n=0 ^ 

Since Yi(x) — xt + 0(t 3 ) is a power series in t with positive polynomial coefficients in x, it 
follows that Y n (x) = xt n + o(xt n ). Hence we have that lim Y^(x) = as a formal power 

iV— >oo 

series in t and the theorem follows. □ 



2.3 The abundant singularities of W(t) 

Our path to proving Theorem 11.11 is now clear. We demonstrate an infinite set of singu- 
larities for W(t) = Q(l, 1) coming from the set of non- conflicting poles of Y n (l). Given 
equation (12. 7p . there are 3 potential sources of singularities for Q(l, 1): 

1. the simple pole at t = 1/3; 

2. singularities from the Y n (l;t); 

3. divergence of the sum. 

We now show that the singularity at t = 1/3 is the dominant singularity and that there 
are an infinite number of other singularities (given by (2)), and that the series is convergent 
elsewhere. 

Lemma 2.3. The generating function Ws(t) has a simple pole at t = 1/3. 

Proof. To show that t = 1/3 is a simple pole, we show that the numerator of equation ( 12. 7h . 
1 — 2tQ(l, 0), is non-zero and absolutely convergent at this point. 

We show in the proof of Corollary 12 . 71 that the radius of convergence of Q(l, 0) is bounded 
below by ~ 0.335, and thus at \t\ = 1/3 the numerator converges absolutely. 

Next we show that the singularity is not removable by proving that 1 — 2tQ(l, 0) is non- 
zero at t = 1/3. Note that K (l; 1/3) = 1 and Y x {\\ 1/3) = 1/2. Using the recurrence (1231) 
for Y n we have 

1 _ 3 1 

y n (i ; i/3) ~ i^oTVs) ~ i^- 2 (i; i/3) 
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which leads to Y n (l; 1/3) = 1/F 2n , i.e. the reciprocal of the 2n th Fibonacci number. Hence 
YrXn+i and so the alternating sum Q(l, 1; 1/3) is convergent. Further, we may write 



( — I)" 1 — -Fkfc-2 

n>0 «>1 

2 10 F i kF i k+2Fik+A 

Hence the sum is bounded strictly between 2/5 and 1/2 and 1 — 2Q(1, 0; l/3)/3 is non-zero, 
and thus t — 1/3 is a simple pole of Q(l, 1). □ 

2.4 The singularities are distinct 

Next, we show that this function possesses an infinite number of singularities — coming from 
simple poles of the Y n . In order to do this we make the substitution t i— > j^-^ which allows 
us to write Y n in closed form. 

Lemma 2.4. Define Y n (q) := (Y„, (l; i^r)) ■ ^/ien 

( q -Y x (g) \ ( l-qZjq) 
1-q 2 J q \ 1-q 2 

1 _ n 2n \ / 1 — n 2n - 2 ' 

1 ^ q^-Y^q)- Vt) (2.9) 



— q 2 J J \ 1 — g 2 

Proof. Substitute i = ^q^2 into equation (12.51) . The roots of the characteristic equation are 
simply q and 1/g and the recurrence can be solved using standard methods. □ 

Since we now have Y n in closed form we can determine some facts about the location of 
its poles. 

Lemma 2.5. Suppose q c is a zero ofY n (q) := Y n (\\ jf^j , and that q c ^ 0. Then 

1. q c is a solution of (q 2n + q~ 2n ) + (q 2 + q~ 2 ) = 4; 

2. q c is not a k th root of unity for any k; and 

3. for all k^n, Y k {q c ) ^ 0. 
Proof. If Y n (q c ) = we must have 

? c (i-»(?c) = (g c 2 -e). 
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Using the equation Yi(q) = {\ + q 2 + y/l — 6q 2 + g 4 j /2q, this can be rearranged to give 

(1 - g 2 ") (l + q 2 c + y/1 - 6q 2 + qt) = 2(g 2 - gf), 

(1 - ? c 2r Vl - 6<Z C 2 + ? c 4 = 2(g 2 - gf ) - (1 - gf)(l + g c 2 ) 

= _(1 _ ql)(\ _|_ g 2n^ 

and so, (1 - ql n f{\ - 6g c 2 + g c 4 ) = (1 - g c 2 ) 2 (l + g 2 ™) 2 . 
We expand and collect terms to arrive at the equation 



which can in turn be reduced to 

+ C 2 ") + + C 2 ) = 4. 

It is clear that q = ±1 are solutions of the above equation. However, substituting g = ±1 
into equation (I2.9p then we find Yi(±l) = ±(1 — i)/2 and 

Y n (l)=nY 1 (l)-(n-l) 
F n (-1) = (-ir +1 (nF^-l) - (n - 1)) . (2.11) 

Hence neither equation can be satisfied and q = ±1 are not a poles of Y n . Now if g c = e l9 , it 
follows that cos(2n#) + cos (26) = 2 so that the only possible zeros on the unit circle are at 
g c = ±1, which we have already excluded. All the zeros lie off the unit circle. 

Next, we prove part (3) of the lemma. First consider the possibility that F n+1 (g c ) = 0. 
If this is true, since the sequence of Y n satisfies a homogeneous three-term recurrence, it 
follows that Y n+ k{q c ) = for all k > 0. It also follows that all Y n _ k (q c ) must be zero. 
However, we can compute the zeros of Y\(q) explicitly and the only possible zero is at g = 0. 
Hence either g c = or Y n±1 (q c ) ^ 0. 

Using the recurrence Y k (q c ) = (j-^^j ^fc-i((Zc) ~ ^fc-2(g c ), and the hypothesis that 

Y n (q c ) = 0, we solve the recurrence for the Y n+k (q c ). This gives 



2A- 



Y n+k (q c ) = y^frj -i-gs— (2.12) 



Since Y n+ i(q c ) ^ by our earlier remark, and since q c is not a root of unity, we have that 
Y n+k (q c ) ^ 0. _ 

By iterating backwards one can similarly show that Y n _ k (q c ) 7^ 0. In fact, one can show 
that Y n+k (q c ) + F n _ fc (g c ) = 0. □ 



Next we show that any such zero is a pole of Q ^1, 0; jf^^j, and hence of Q ^1, 1; — 



7. 
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Lemma 2.6. Let q = q c ^ be a zero ofY n . The function Q ^1,0; j+pj has a pole at the 
Q = Qc- 

Proof. First, recall that as = q k + 0(q k+1 ), this is a convergent expression (as a formal 
power series), and thus it makes sense to consider this series in this arrangement. We contend 
that the only poles come from the poles given by zeroes of Y n , for all n. Aside from these 
poles and the square root singularity given by q = =t-/g, the series is convergent. 
Write YfcYfc+i as 

Y t Y M = (1 - f) (1 - Y iq) ( q , HY ._ l q) _ Yi+q - , M( , ri _ 1 1) _ yl+g ) ■ (2.13) 

Suppose q is such that the denominator is not zero, and such that it does not lie on the 
branch cut of Y\. 

We know that \q c \ ^ 1, so we only need to consider \q\ > 1 and \q\ < 1. The above 
expression for YkYk+i is obtained by substituting t \— > and so is necessarily invariant 
under q t— > 1/q. Hence we only need to consider \q\ > 1. In this case the denominators in 
the right-hand side of the expression are dominated by the q k term and so approach as 
k — > oo. Since the sum is alternating it converges. 

Now, the denominator is annihilated in this case precisely when q is a zero of either Y^{q) 
or Y k+ i(q). Suppose q c is a zero of Y k (q). Then, by Lemma |275| Y n (q c ) ^ if n ^ k. Thus, 
the only two terms for which q c is singular are Y k _\f k and Yfelfc+i. If we sum these two 
terms we have (— l) h Yk{Yk + i — ife-i). By our earlier comments, Yk-i(q c ) = —Yk+i(q c ) ^ 0, 
and thus this singularity does not cancel. □ 

2.5 Proof of Theorem 11.11 

We now have all the components in place to prove the main result. 

Proof of Theorem ] 1 . 11 The function Q ^1, 1; jf^j has a set of poles given by the zeroes of 

the Y n , by the preceding lemma. By Lemma [231 this is an infinite set. Thus, Q ^1, 1, is 
not holonomic. For a multivariate series to be holonomic, any of its algebraic specializations 
must be holonomic, and as Q ^1, 1, jf^ij is an algebraic specialization of both Q(x,y) and 
(5(1, 1), neither of these two functions are holonomic either. □ 

Corollary 2.7. The number of random walks of length n with step set {NE, SE, NW} confined 
to the quarter plane is asymptotic to 

c n ~ aT + 0(8 n/2 ) (2.14) 

where a is a constant given by 

(—l) n 

ft = 1 " 2 EF7 — = 0.1731788836- •• (2.15) 

n>0 F ^F2n+2 
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Proof. We need to show that the dominant singularity is the simple pole at t = 1/3. Let us 
rewrite equation (12 .ip at x = y = 1: 

(l-3f)Q(l,l) = l-2«Q(l,0) or Q(l, 1) = 1 ° } . (2.16) 

Hence there are 2 sources of singularities for Q(l, 1) — the simple pole at i = 1/3 and the 
singularities of Q(l, 0). The residue at the simple pole may be computed quite directly (see 
the proof of Lemma 12. 31) . 

We now bound the singularities of (5(1,0) away from \t\ = 1/3. The series Q(1,0) 
enumerates random walks confined to the quarter plane which end on the line y = 0. Consider 
a family of random walks with the same step-set but ending on the line y = but no longer 
confined to the quarter plane, rather just confined to the half-plane y > 0. The generating 
function of these walks obeys the functional equation 

P(y) = l + t[2y + -) P(y) - t-P(O). (2.17) 

V yj y 

The generating function of these walks can be computed using the kernel method to give: 



Hence the number of these walks grows as 0(8 n / 2 ). Consequently the walks enumerated 
by Q(1,0) cannot grow faster than 0(8 n//2 ), and so the radius of convergence of (5(1,0) is 
bounded below by 8 -1 / 2 . The result follows. □ 

3 The non-holonomy of Qt{x, y\ t) 

This step set is not symmetric across the line x = y, and thus as we lose x <-> y symmetry 
in the complete generating function, we expect a more complicated scenario. Nonetheless, 
we can recycle a good number of the results from the previous case. 

3.1 Defining two sets of iterates 

We begin by rewriting the fundamental equation in kernel form: 

K(x, y)Q{x, y) = xy- txQ{x, 0) - ty 2 Q{0, y), (3.1) 

with kernel K(x,y) = xy — t(xy 2 + x 2 + y 2 ). We will require zeros of the kernel for both 
fixed x and fixed y, which we denote respectively with X and Y. This gives, 

X±i(v) = Y t i lT V(l-ty) 2 -t 2 ) , and (3.2) 
Y±i(x) = (l =F y/1 - U\x + 1)) • (3.3) 
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Using similar argument to those in Section [2j one can show that 

X ± (Y T (x))=x and Y ± (X T (y)) = y (3.4) 

and that 

1 i-l (3.5) 



X+i(y) X-ifo) ty 

1 + ^ = ^. P.6) 



We can use these to find nice expressions for compositions of X's and Vs. Define B n (y) and 
A n (x) by 

Bo(y) = y 

B 2k+1 (y) = X 1 (B 2k (y)) (3.7) 
B 2k+2 (y)=Y 1 (B 2k+1 (y)) 

and 

Aq(x) = X 

A 2k+1 (x)=Y 1 (A 2k (x)) (3.8) 
A 2k+2 (x) = X 1 (A 2k +i(x)). 

As before, we can compose these relations and determine the recurrences 

11 1 1 



B n+ i(y) tB n (y) 
1 1 1 



A n -\-\(x) A n _i(x) tA n {x) 



(3.9) 



Note that A k satisfies the same recurrence and initial conditions as Y k in Lemma 12.11 of 
Section [2] 

3.2 Expressing the series in terms of the iterates 

Next, we use the relations such as K(A 2k ^i(x), A 2k (x)) = and K(A 2k (x), A 2k+ i(x)) = to 
generate an infinite set of equations. For example, if we set y = Ai(x) = Y(x), the functional 
equation Eq. (13.11) reduces to 

Y x (x) (Y 1 { X y 2 



Q(x,0) = -^--l-^) Q&Y^x)), 



and when x = B\(y) = X\(y) 



X 1 (y) fX 1 (y) 



2 



Q(0,y) = - Q(Xi(y),0). 

yt \ y 



n 



Putting y = Yi(x) into the above expression for Q(0,y) gives 



Xi(Yi(x)) fX 1 (Y 1 (x)) 



2 



Rewriting everything in terms of A n (x) and substituting one equation into the other we get 

= m^mi _ c 4441 _ M&ay Q( ^ Wi0) 



Ao(x) 2 t \A Q (x) J \ Ai{x) t \A 1 {x) 



A (a;) 2 t VA)0 

Iterating this procedure leads to the following theorem 

Theorem 3.1. The generating function for the walks that end on either axis and the counting 
generating function are given by the following expressions: 



Q(x,0) = -^-T(-l) n A n (x)A n+1 (x), 
art L — ' 

n>0 

Q(°> y } = \ 12^ nB M B ^(y^ ( 3 - 10 ) 

" n>0 

0(1, i) = (1 - ^(-i) n A n (i)A n+1 (i) - ^(-i) n B n (i)fi n+1 (i) J . 

\ n>0 n>0 / 



3.3 A source of singularities 

We now show that the generating function has an infinite number of poles, and that these 
are given by the A n , and the B n . We consider the series under the same g-transformation 
that we used when we considered the other step set. We set A n (q) = A n (l; jf^)" 1 and 
B n (q) = B n (l; T ^y\ 

Lemma 3.2. The function A n (q) is identical to Y n (q). 

Proof. When x — 1, the kernels for the functional equations for both step sets are the same. 
This implies that when x = 1, the solutions to the respective kernels are the same; That is, 
Y + is the same for both step sets and Y~_ is the same for both step sets. 
Consider the recurrence satisfied by the A n . Iterating Eq. (13.91) gives 

t/x fl-q 2n \ 1 -r , x fl-q 2n - 2 \ 1 , , 

= { T^qT) ^ Mq) ~ I^VJ ^' (3 - U) 

and thus, A n (q) here, satisfies the same recurrence as Y n (q) given in Eq. ( 12.121) . and also the 
same initial condition. Thus we conclude, A n (q) = Y n (q). □ 
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Corollary 3.3. Let q = q c 7^ be a zero of A n (q). Then, 



1. q c is not a zero for A n (q) whenever n ^ k; 

2. The function Q (1,0; ) has a pole at q = q c . 



Proof. This follows from the previous lemma and Lemmas 12.51 and 12.61 



□ 



As in the case of step set S, we have that tQ (1,0) = ^2(—l) n A n (l;t)A n (l;t) has an 
infinite number of singularities and so in not holonomic. Thus consequently Q(x,y), is not 
holonomic (since one is a specialisation of the other). Thus, already we have proven the first 
part of Theorem 11.21 

In order to prove that the counting function W(t) = Q(l, 1) is also not holonomic, we 
must show that the singularities of Q(l, 0) are singularities of W(t) and so are not cancelled 
out by the singularities of Q(0, 1). This is what we do next. 

3.4 The singularities do not cancel 

We remark that we can deduce an expression for B n B n+ i similar in spirit to that of Eq. (I2.13p . 
Thus, for all points which are not zeroes of some B n (q), the sum ^2 k>0 (—l) n BkBk+i con- 
verges. 

Now, the singularities for Q(1,0) and Q(0,1) have different sources, and thus a priori 
we expect that they will not be the same. We have not been able to do this. However, it 
is sufficient (and easier) to show that there is an inifite subset of singularities that do not 
cancel. In particular we consider the singularities of the A n that lie on the imaginary axis. 

For even n we will show that A n has poles along the imaginary axis. Now B n has any 
pole along the imaginary axis, and thus Q(l,l) has an infinite number of singularities along 
the imaginary axis. We prove this technical lemma now, and with that our work is done. 

Lemma 3.4. The functions A n (q) and B n (q) satisfy the following recurrences: 




(3.12) 



.4 



Aq 



A 



1 



■ 2 + l + vV - 6q 2 + 1 




B 



2q 



B, 



1. 



q 2 - q + 1 + vV - 2q 3 - q 2 - 2q + 1 



o — 



Thus, it follows 



(1) For n >1, A2n has at least two zeros along the imaginary axis; 
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(2) For any n, and any q = ir along the imaginary axis the value of B n is non-zero, except 
at q = 0. 

Proof. Equation (I3.12p is a direct consequence of Equation (I3.9p . (1) From the proof of 
Lemma 12.51 in Section [2] we see that the zeros of Y n are given by the zeros of the following 

equation 

(1 - q 2n ) y/1 - 6q 2 + g 4 + (1 - q 2 ){l + q 2n ) = 0. (3.13) 

By Lemma 13.21 we know that the zeros of A n satisfy the same equation. 
Now substitute q = ir and n = 2m into this equation to obtain 



(l-r 4m )Vl + 6r 2 + r 4 + (1 + r 2 )(l + r 4m ) = 0. (3.14) 

Denote the left-hand side of this equation by f(r), which is a continuous function for real r. 
At r = ±1 we find that /(±1) =4. At r = ±2 we find that 



f(±2) = (V41 + 5) + (5 - v / 41)2 4m (3.15) 

which is negative for m > 1. Hence f(r) has at least two zeros, one in the interval [—2, —1] 
and the other in the interval [1,2]. So A2 n has zeros on the imaginary axis. 
(2) First, we solve the recurrence and give an explicit form for B n in terms of B\, denoted 
here by /3: 

- _ q^(l - 2 q + (3q 2 - (3q) + q n (q 2 + q) + (g 3 - 2q 2 + q(3 - (3q 2 ) 

(-l+q) 2 (q+l)qn 

We can now repeat the method used in the proof of Lemma 12.51 to show that the zeros 
of B n satisfy the following equation: 

= q 2 (q 2 -3q+l) {q in + l) - q(q + l) 2 (q 2 - 3q + 1) (q 3n + q n ) 

- ( g 6 _ 4g 5 + 3g 4 + Qq 3 + 3q 2 - Aq + l)q 2n . (3.17) 

Note that this equation is symmetric in q <-> 1/q (as we would expect, since we introduced 
the q variable by substituting 1 1— > l/(q+ 1/q))- 

We wish to show that the solutions of this equation do not lie on the imaginary axis and 
so the singularities of A n which lie on the imaginary axis cannot be cancelled by singularities 
of B n . Hence we substitute q = ir with rGK: 

= r 2 (r 2 + 3ir - 1) (r 4n + l) - ir(r - i) 2 (r 2 + 3ir - 1) {(irf n + (ir) n ) 

+ (r 6 + 4ir 5 - 3r 4 + 6ir 3 + 3r 2 + Air - l)r 2n (-l) n . (3.18) 

Let u(r) denote the right-hand side of this expression. Note that 2/(0) = for all n and also 
z/(±l) = ±81, so that r = is always a zero and r = ±1 are never zeros. 

We now analyse the real zeros of this expression for different values of n mod 4. In each 
case we show that there are no real values for r which annihilate v (except r = 0). We do 
this by considering the real and imaginary parts of v separately. 
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Case 1: n = 2 (mod 4) 

In the case that n = 2 (mod 4), the imaginary part of v is 

= 3r 4n+3 + r(r 4 + 4r 2 + l)r 3n + 2r(2r 4 + 3r 2 + 2)r 2n + r(r 4 + 4r 2 + l)r n + 3r 3 . (3.19) 

This expression is strictly positive for r > 0. Since it contains all odd powers of r, it is 
strictly negative for r < 0. It can only be zero when r = 0. 

Case 2: n = (mod 4) 

In this case, v has real part equal to 



implies that r 2n + r™ + r 2 + r~ 2 + r _n + r _2n = 2. This equation is symmetric in r <-> 1/r 
and r <-> — r (since n is even). If r > 1, the left-hand side is strictly larger than 2, and hence 
has no solution. By symmetries, there can be no real solution. 

So the only points that annihilate the real part of v are r = and r = ±1. As noted 
above r = is a zero of is, but r = ±1 are not. 

Case 3: n odd 

To show that is has no real zeros (aside from r = 0) when n is odd, we consider a linear 
combination of the real and imaginary parts of u that eliminates the leading powers of r 
which in turn simplifies the analysis. 

In particular consider (r — l/r)Q(u) — 3dt(is); this eliminates the (r 2 + 3ir — l)(r 4n + 1) 
term from equation (13.181) . If there exists a real zero of v then it must also annihilate this 
linear combination. We will show that no such value of r exists, and so there are no real 
zeros of v (excepting r = 0). 

The linear combination simplifies to 



(r - l/r)Q(u) - 38l(u) 

= r 2k+2 (2 + Ur 2 + 2r 4 )(r 4fe+2 - l)(-l) fc - (r 2 - l)(r 4 + 12r 2 + l)r 4fc+2 . (3.22) 



Now divide both sides by r 2 — 1 to obtain (we know that r ^ ±1 are not zeros of is): 



= r 2 (r 2 - 1) (r 4n + r 3n + (r 2 - 2 + r" 2 )r 2n + r n + l) . 



(3.20) 



Suppose r 7^ 0, ±1. Then 



r 2 (r 4n + r 3n + r n + r + (r 2 - 1)V) = 



(3.21) 



r 



1 



) 



(2 + 14r 2 + 2r 4 ) 



(_!)* _ ( r 4 + 12r 2 + i) r 4fc+2. (3 23) 



r 



2_! 
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Figure 2: Stretching the walk to find a directed path in a strip 



When k is odd this expression is a polynomial in r 2 with negative coefficients, while when k 
is even this is a polynomial in r 2 with positive coefficients. Hence this is non-zero except at 
r = 0. 

Thus, as no B n has a zero on the imaginary axis, (2) holds. □ 

It seems clear that Q(0, 1; t) is also non-holonomic. Indeed, the walks are very similar to 
the symmetric case. However, we were unable to construct a direct argument to show that 
it too is non-holonomic. 

4 Conclusions 

4.1 How robust is this method at detecting non-holonomy? 

This is a natural question, and it speaks to the ability of this approach to be applied to other 
problems. In all the cases known to give a holonomic generating function, the group of kernel 
iterates is finite (of order 8 or less). If one attempts to repeat the proofs of Theorems 12.21 
and 13 .11 then one obtains a finite set equations. It is not clear that solution of these equations 
has singularities appearing in the same way as occurred here. 

4.2 Combinatorial intuitions of non-holonomy 

More generally, we are interested in developing an intuition about when precisely a combi- 
natorial object will have a non-holonomic generating function. If we decompose these walks 
according to their NE or N steps, a source of singularities becomes apparent. Unfortunately, 
since we are unable to prove that the singularities do not cancel, the details we present in this 
section do not constitute a proof of non-holonomy, however they are nonetheless instructive 
for understanding a potential source of the singularities. 

Any walk from either of these two families can be decomposed into a sequence of directed 
paths in a strip of unrestricted length where after each N E/ N step the size of the height bound 
increases. Figure [2] for an example of such a decomposition for step set S. Thus, a potential 
strategy groups lattice walks that end on the line x + y = k, and decomposes any such walk 
into a triple: A walk that ends on the line x + y = k — 2, a NE step, and a directed path in 
a strip of height k. 
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To describe the generating function, we use Example 11 of [5] who present a generating 
function for paths of length n in strip of height k of walks that begin at a given height, and 
end at a given height. 

Define D k (y; t) as the generating function for the subset of walks ending on x+y = k where 
y marks the final height of the walk. We can easily translate the above decomposition into a 
functional equation for the generating function. If we then make the substitution t 
the expression simplifies remarkably into the following recurrence for D k (y) = D k (y 



1+9 
_2_^ 



2 • 



' 1+q 2 



q 3 D k . 2 (q)(y k+2 + 1) - qy 2 D k - 2 (y)(q k+2 + 1) 
In fact, for our purposes it suffices to consider: 

Dk{1) = (<f+ 2 + l)(q - l) 2 • (42) 

From this formula, and from computations for various values of k, D k (l) is a rational function 
in q, and it seems clear that the set (taken over all k) of poles of D k (l) is dense in the 
unit circle. Were this so, we would apply the following theorem to the generating function 
Qs(s,s; jf^z) = £)-Dfc(l; jf^)s k , and thus conclude the non-holonomy of Qs(x,y;t). 

Theorem 4.1. Let f(x;t) = 'Yl, n c n {x)t n be a holonomic power series in C(x)[[t]]. with 
rational coefficients in x. For n > let S n be the set of poles of c n {y), and let S = [jS n . 
Then S has only a finite number of accumulation points. 

Again, the principle difficulty is showing that the singularities do not cancel; that solu- 
tions to q k+2 + 1 are indeed poles of D k . 

This approach was pioneered by Guttman and Enting [7J , and has been fruitful for several 
different models [HI [12]. Unfortunately their arguments do not appear to work here. 



4.3 Related walks 

We expect the following walks to have non-holonomic generating functions because the 
groups of their kernel iterates are infinite. It is even likely this can be proved in the same 
manner as Theorems 11.11 and II .21 
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